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Motivated by several experimental efforts to understand spin diffusion and transport in ultracold 
fermionic gases, we study the spin dynamics of initially spin-polarized ensembles of harmonically 
trapped non-interacting spin-1/2 fermionic atoms, subjected to a magnetic field gradient. We obtain 
simple analytic expressions for spin observables in the presence of both constant and linear magnetic 
field gradients, with and without a spin-echo pulse, and at zero and finite temperatures. The 
analysis shows the relevance of spin-motional coupling in the non-interacting regime where the 
demagnetization decay rate at short times can be faster than the experimentally measured rates 
in the strongly interacting regime under similar trapping conditions. Our calculations also show 
that particle motion limits the ability of a spin-echo pulse to remove the effect of magnetic field 
inhomogeneity, and that a spin-echo pulse can instead lead to an increased decay of magnetization 
at times comparable to the trapping period. 
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I. INTRODUCTION 

Understanding spin transport in quantum systems is 
central to many fields of physics and can have impor¬ 
tant applications for the development of quantum tech¬ 
nologies. Recently there has been much experimental 
progress studying spin diffusion m and spin segrega¬ 
tion dynamics (time-dependent separation of the spatial 
distributions of the spin components) [TUHTS] in trapped 
atomic gases. A typical experimental protocol consists 
of preparing a transversely spin-polarized gas and then 
observing the spin relaxation dynamics under the influ¬ 
ence of a magnetic field gradient. For instance. Ref. [5] 
measured the demagnetization timescale in a constant 
magnetic field gradient and determined the spin diffu¬ 
sion constant. Similarly, Ref. m measured the segrega¬ 
tion of the spin populations in a magnetic field gradient 
with linear spatial dependence. Although the goal of the 
experiments is to understand the many-body interact¬ 
ing spin dynamics, it is important to have a clear under¬ 
standing of the non-interacting physics and how the spin- 
motion coupling alone affects the spin demagnetization. 
A thorough analysis of the non-interacting system will 
help provide the foundation necessary for understanding 
the complex spin dynamics induced by the interplay be¬ 
tween interactions and motional effects. 

Here, we provide analytic expressions for the demag¬ 
netization exhibited by a harmonically trapped and non¬ 
interacting spin-1/2 Fermi gas at both zero and finite 
temperatures. We study the spin dynamics with and 
without a spin-echo pulse and in the presence of mag¬ 
netic field gradients with constant or linear spatial de¬ 
pendence. The spin dynamics is oscillatory, and depends 
on details of the differential motion of the spin compo¬ 
nents. In the case of a constant gradient (magnetic field 
that varies linearly with position), the atoms oscillate at 
the trap frequency but in opposite directions along the 
gradient depending on their spin projection. We show 


that this periodic motion gives rise to a fast demagneti¬ 
zation but not to a net spin rotation. In the case of a 
linear gradient (magnetic field that varies quadratically 
with position), the spin components breathe at different 
rates and with different magnitudes, the dynamics is not 
periodic with the trap frequency, and the spin dynamics 
involves both a demagnetization and a net rotation of the 
collective spin. 


In the non-interacting limit our analysis reveals that 
the transverse magnetization decays with a rate that 
grows with increasing particle number and temperature. 
The observed fast demagnetization rate at short times 
reflects the fact that samples with a large number of 
fermions occupy high harmonic oscillator modes that are 
widely spread across the trap and experience strong gra¬ 
dients. We also find that the spin-motion coupling can¬ 
not be removed by a spin-echo pulse, and such a pulse 
can instead lead to an enhanced decay rate of the mag¬ 
netization of the gas. In unitary Fermi gas experiments 
this fast motion-induced demagnetization is suppressed 
at short times by interactions which instead lead to a 
diffusive decay of the magnetization. 


This paper is organized as follows: In Sec.|TT]we present 
and discuss the spectroscopic protocol under considera¬ 
tion. In Sec. in and Sec. |IV] we present the constant 
and linear gradient cases, respectively. In each of these 
sections we derive analytic expressions for the single par¬ 
ticle dynamics and then formulate the spin dynamics of 
the many-particle system at finite temperature in terms 
of these expressions. We derive expressions for the dy¬ 
namics both in the presence and absence of a spin-echo 
pulse and discuss experimental considerations. Finally, 
in SecjVj we finish with an outlook and conclusions. 
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II. SPECTROSCOPIC PROTOCOL 

We begin by considering a single spin-1/2 particle of 
mass m confined in a one-dimensional harmonic trap with 
trapping frequency uj - results for ensembles of atoms are 
later calculated from sums over single particle dynamics. 
The spin dynamics is probed using Ramsey spectroscopy: 
at t = 0 the particle is prepared in an eigenstate of the 
harmonic trap and with spin j/). The spin of the particle 
is then rotated about the X-axis by applying a resonant 
pulse with area 9i. Next, a position-dependent magnetic 
field Bz{x) pointing along Z is suddenly turned on, and 
the particle then evolves freely in the presence of the mag¬ 
netic field for a dark time t, after which spin observables 
are measured. 

Here, we focus on the case of Ramsey spectroscopy 
with an initial pulse area 9i = tt 12 that rotates the initial 
state to point along Y (transverse direction). We also 
assume that the pulse has zero detuning from the atomic 
transition (5 = 0). However, the results that follow are 
easily generalized to arbitrary 9i and finite detuning. 

During the dark time the particle feels a potential 

V{^) = V{x)+V{y,z) 

V{x) = + AfiBz{x)a^, (1) 

where A/r is the differential magnetic moment between 
the two spin states and is the usual Pauli operator. 
Along the y and z directions the potential is assumed 
to be spin-independent. The dark time evolution of the 
state of the particle can be written as 

|^'(x,t)) = t) + V^t(x,t)| i). (2) 

Due to the separability and the spin independence of the 
V{y, z) potential, if the system is prepared in an eigen¬ 
state of y and z, it will remain in that eigenstate and thus 
the dynamics is effectively one-dimensional. Because of 
this, we will restrict our analysis to the x dimension. 

We focus on the collective observable + iS ^, 

where 5“ are the spin angular momentum operators. The 
expectation value of S~^ takes the form of an integrated 
density 

(* 5 +) = ^ J dx'tlj^x,t)ipi{x,t) 

^ (3) 

Here the Ramsey fringe contrast C{t) = |A(t)|/|A(0)|, is 
related to the overlap of the ip^{x,t) and ipi{x,t) wave- 
functions. The decay of the contrast is a measure of de¬ 
magnetization. The frequency shift Av{t) measures the 
dynamical phase difference between the spins and gives 
rise in this case to a net motionally-induced precession 
of the total magnetization. Throughout we set h = 1. 


III. LINEAR MAGNETIC FIELD 

First we consider the case of a magnetic field with lin¬ 
ear spatial dependence Bz{x) = Bx/Ay, where H is a 
constant with units of energy/length. Adding a linear 
potential to a harmonic potential results in a new har¬ 
monic potential with the same frequency but shifted in 
opposite directions for each of the two spin states. The 
potential is 

V{x) = ^muj^{x + xq)'^ — ^moj'^XQ . (4) 

Here, an = \/llmu: is the oscillator length and xq = 

is the displacement of the oscillator resulting from 
the magnetic field. 

It is well known that the displaced ground state of 
a harmonic oscillator evolves with a probability dis¬ 
tribution of constant shape but oscillating centroid, 
11 / 0 ( 2 :, t)P oc [16]. A similar result can 

be derived for any eigenstate n. Namely, given a solution 
4>[x, t) to the time-dependent Schrodinger equation, 

^{x, t) := b (zoe"“*) (j){x, t), (5) 

is also a solution to the time-dependent Schrodinger 
equation. The displacement operator is given by D(w) = 
exp (w(b — w*d) and Zq is any complex number, cor¬ 
responding to an initial displacement in the position- 
momentum phase space mm- This result allows us to 
calculate dynamics analytically and write the time evo¬ 
lution of an initial eigenstate as (See Appendix A): 
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—if n+ ^ ^ cos(u;i) sin{ujt)—sin{ujt )) 

(X.t) =e 

„ 31 + aXp (1 - COS (o^t)) ^ ^ 

where a is the eigenvalue of . In Fig. [^a), we plot 
the time evolution of the spin density in the Z-direction, 
{S^{x,t)) = i (|V)^(a:,t)p - |i/;(a;,t)p), for a single par¬ 
ticle initially in the ground state of the harmonic oscil¬ 
lator, calculated using Eq. The | f) and | 1) densities 
maintain their same shape but oscillate about their re¬ 
spective equilibrium positions at the trap frequency. 

A. Spin dynamics 

We now use the above results to determine the spin 
dynamics for a single particle. Assuming that the particle 
begins in the harmonic oscillator mode at time t = 0, 
we find 

0b = \ exp ^-2^ (1 - cos(wt))^ 
xL„ ^4^(1 - cos{ujt))\ , 


( 7 ) 
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FIG. 1. Dynamics of spin observables in a linear 
magnetic field with xo = an- (a) {S^{x,t)) = 

I shows the spin up/down densities 

oscillating in their traps centered at a; = =F®o. (b) Magnitude 
of {S'^ {x,t)), which measures the magnetization at point x. 
The magnetization decays when the spin up/down densities 
are separated. 


where L„ is the n**' Laguerre polynomial. 

We see that the contrast is suppressed exponentially 
with increasing xq. This is caused by the exponential 
suppression of the overlap integral between the spatial 
wavefunctions of the f and | spin states. As they move 
apart under the influence of their different potentials, the 
overlap decreases as a consequence of the wavepackets’ 
Gaussian spatial localization, see Fig. int)- The contrast 
is periodic in time with period 27r/a;, due to the periodic 
motion of the two spin wavefunctions in the trap. We also 
see that there is no frequency shift for the case of a lin¬ 
ear magnetic field and the direction of the magnetization 
remains along Y. The two spin states, while displaced 
in position, experience harmonic potentials of the same 
frequency, so they have the same dynamical phase. 

We now consider the dynamics of an ensemble of non¬ 
interacting fermions with this same protocol. An ensem¬ 
ble at zero temperature, initially spin-polarized, forms a 
Fermi-degenerate gas with all oscillator modes filled from 
n = 0 to n = N—1, with N the number of particles. Here, 
we find 


N-l 


{S 


+ \f.d. 




n—O 


i —2—^(1—cos(cjt)) 
= -e 


N-l 


E Ln ( 4^(1 - cos(w<)) 

n=0 k Off 


2® 


i —2^-(l—cos(wt)) ^ 


_P 
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^iV-i 4^(l-cos(a;t)) , (8) 


where L'^(x) is the associated Laguerre polynomial and 
we have used the addition formula for Laguerre polyno¬ 
mials Em =0 ^n-m iv) = {x + v) ■ The 

contrast decay per particle at short times is given by 


C{t)/N = 1- Ep^ojt^ + O (f^) , 


( 9 ) 




with Ep = ksTp = Nuj the Fermi energy of the gas, kp is 
the Boltzmann constant, and Tp the Fermi temperature. 
Hence, the decay of the contrast per particle increases 
extensively as the number of particles increases. 

For an ensemble at non-zero temperature, the associ¬ 
ated Fermi-Dirac sums cannot be evaluated exactly, but 
can be approximated by energy integrals with a constant 
density of states to yield a formula for the short-time 
contrast decay per particle: 

C{t)/N = l- 2E^ojf + O (t^) 

Qh 

TT^ + 3 ln(e^'^“ - 1) + 6Li2 , 

( 10 ) 

where P = l/kpT, T the temperature and Lig {x) is the 
polylogarithm of order s. 

In the limit of high temperature, T/Tp ^ 1, we can ap¬ 
proximate the distribution function of the fermions with 
a Maxwell-Boltzmann distribution. For this ensemble the 
spin dynamics is given by 


E 


6NioP^ 


n—0 

1 Xq (1 — cos(a;t)) 

2 ^ tanh(/3a;/2) ’ 


( 11 ) 


where Z = e is the partition function. 

At short times, we find the contrast per particle de¬ 
cays asCit)/N=l-^ tanhE/ 2 ) ^"^" + ^ ’ ^hich 

agrees with Eq. ( |Io|) when we identify the mean energy 
1 tanh“^(/3tLi/2) = E/co. At T/Tp 1, where this anal¬ 
ysis is valid, E « and so the contrast decay rate 

increases linearly with temperature. 


B. Dynamics for a spin-echo sequence 

The spin-echo consists of an additional collective tt ro¬ 
tation about X added at time t/2. This pulse swaps the 
states of the spin up and spin down components, with 
the goal of removing spin-dependent single-particle inho¬ 
mogeneities. It is natural to wonder whether a spin-echo 
pulse will remove the effects of a magnetic field gradient, 
which is effectively a spatially-inhomogeneous detuning, 
when the particles themselves move in the trap. In the 
presence of a spin-echo pulse the evolution is given by: 

(SysE = iexp(^- 16 isi„‘(^)) 
xL„(32;|si„‘(:^)j. (12) 

The spin dynamics with a spin-echo pulse is periodic with 
period 47r/(jj, twice the period of a Ramsey sequence with¬ 
out a spin-echo pulse. The short-time expansion of the 




























contrast decay is 
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CsE (^) = 1 - § (^ + ^) ^ (^®) • (13) 

Note that the spin-echo removes the dominant O (t^) 
contribution to the contrast decay at times much shorter 
than the motional period. This is consistent with the 
expectation that spin-echo removes the effect of inhomo¬ 
geneities when motional effects can be neglected [H]. 
Beyond the short time limit the spin-echo pulse is not 
beneficial. 

The behavior of the spin-echo signal at times compa¬ 
rable to the trap period can be visualized in phase space. 
To illustrate the dynamics, we rescale the phase space 
coordinates to x/xq and p/po, where po = I/oq, see Fig. 

g . In this coordinate system, the expectation values 
and (p) for the | f) state and | 1) states initially 
move along circular trajectories centered at ^xq, reflect¬ 
ing their initial displacement from the center of their re¬ 
spective harmonic traps due to the magnetic field. The 
spin-echo tt pulse switches the | f) and | j,) states and 
while doing so it generally enhances the net displacement 
of the spin states from their motional centers. After the 
pulse the dynamics corresponds to circular trajectories 
with a new phase-space radius. For example, as shown in 
Fig. i if the echo is applied at tt/w (dark time t = 21 ^ ju))^ 
the displacements following the spin-echo pulse are twice 
the displacements before the pulse. 

For dark times less than half the motional period, 
the spin-echo improves the contrast, consistent with the 
short-time analysis. In the phase space picture this 
means that the displacement between | f) and | j,) with 
spin-echo is smaller than without spin-echo. At a dark 
time of exactly half the motional period, the contrast 
decay with and without spin-echo is identical - this is 
illustrated in Fig. In this case the phase space dis¬ 
placement without spin-echo is purely along x and the 
displacement with spin-echo is purely along p. For longer 
dark times, spin-echo increases dephasing, since it pro¬ 
duces larger phase space displacements. The dephasing is 
maximal at dark times which are odd-integer multiples of 
the motional period. When the spin-echo pulse is applied 
at an integer multiple of the motional period, i.e. dark 
times of integer multiples of twice the motional period 
t = A-nk/uj [k an integer), both spin states have returned 
to the phase-space origin and the motional dynamics is 
unaffected by the spin-echo pulse. This explains the pe¬ 
riodicity of the spin-echo signal with twice the motional 
period. We stress that the ability to decouple the spin- 
echo effect from motion is due to the two spin states shar¬ 
ing a common motional frequency. Although spin-echo is 
typically used to eliminate dephasing due to single par¬ 
ticle inhomogeneities, in the presence of spin-motional 
coupling the particle motion during the dark time can 
worsen the resulting dephasing when a spin-echo pulse is 
applied. 

We can calculate the dynamics of thermal ensembles 
with a spin-echo pulse. At zero temperature, the exact 



FIG. 2. Phase space plots of {x) and (p) for the | 4) state 
(blue, solid) and | 4) state (red, dashed) for a spin-echo se¬ 
quence in a linear magnetic field. Here po = 1/fflo, and the 
total time is 2tx jbji. The two spin states evolve along circular 
trajectories centered at =F®o. (a) At half the total dark time, 
tt/cu, the spin-echo tt pulse swaps the spin states. The spin 
wavefunctions are thus displaced by 2a:o from the centers of 
their respective harmonic traps. For the second half of the 
dark time they evolve in circular trajectories twice as large, 
centered at ^2xo, resulting in a final configuration shown in 
(b). The spin-echo pulse leads to greater dephasing of the 
spin observables due to larger final phase-space displacement. 


result is 


('S'+)sE = o exp ( sin"* ( — 




X A]v-i 


32^sin^ ( — 


(14) 


For an ensemble at arbitrary temperature the short time 
contrast decay is given by 


CsE (t) /N = l- + O (t®) 




E 


6A^w/32 


+ 31n(e^^‘^-l)-b6Li2 


1 


^ _ q^Nuj 


(15) 


and at high temperature we can approximate the ensem¬ 
ble by a Maxwell-Boltzmann distribution to obtain 


(^+)rE^' = ^exp 


-16 


^gsin^(f) 

tanh(/3w/2) 


(16) 


The contrast at zero temperature and high temperature 
(T 3> Tp) are shown with and without spin-echo as a 
function of dark time in Fig. 


C. General Considerations 

In general terms the contrast dynamics in the non¬ 
interacting regime is characterize by oscillatory behavior 
at the trap period arising from spin-motion coupling. At 
short times the contrast decays as (7(1:) « 1 — (|:) and 
in the presence of an echo-pulse as (7®^(1) « 1 — 
where we identify r and r®® as the demagnetization 
timescales without and with spin echo. Using Eqs. 
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and 15 the demagnetization timescales are given by: 


Tiin = — {2i^E) 
Xo 






Min 


E 


Xt\ ^ 


,3 -^0 




- 1/4 


6iVa;/32 


TT^ +31n(e^^‘^ - 1) + 6Li2 


1 

^ _ q^Nlo 


(17) 


FIG. 3. Contrast decay for the case of a linear magnetic field 
at times comparable to the trap period, for zero temperature 
and high temperature {T ^ Tf), with and without spin-echo 
(SE). The quantum degenerate case uses = 19 and all 
cases use xq = 0.24a_f/, consistent with Ref. [2]. T = lOTp 
in the high temperature case. Higher temperature leads to 
faster contrast decay since the particles have a larger average 
energy. The spin-echo removes the second-order contribution 
to the contrast decay, resulting in slower decay at short times. 
However, at long times the contrast decay is larger when a 
spin-echo pulse is applied. The period of the dynamics is also 
4Tr/uj with a spin-echo pulse instead of 27r/a;. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 


T/Tf 

FIG. 4. Linear magnetic field: (a) demagnetization time vs. 
magnetic field gradient, (b) Demagnetization time vs. tem¬ 
perature where the Fermi temperature Tp = Nu). Using pa¬ 
rameters from Ref. [^, these timescales are faster than those 
obtained in the experiment, suggesting that strong interac¬ 
tions suppress the role of single particle dynamics. 


The short time demagnetization rate can be significantly 
faster than demagnetization rates measured in recent ex¬ 
periments performed at unitarity m, where the magne¬ 
tization decay is manly diffusive with a short time scal¬ 
ing given by C'^(t) « 1 — (;^)^. The different time 
dependence translates into a different dependence on 
the magnetic field gradient: while in the non-interacting 
regime the demagnetization rate scales as 1/t ^ B and 
- VB, in the interaction dominated regime it 
scales as I/r ^ [IHS]. 

The demagnetization rate in the non-interacting sys¬ 
tem increases with temperature at fixed particle number 
with and without echo. The increasing decay rate in the 
non-interacting limit simply reflects the increase in the 
mean energy per particle as the gas gets hotter. This 
behavior persists in the presence of interactions as ob¬ 
served in current experiments given that both the typical 
velocity and the mean free path (parameters that deter¬ 
mine the spin diffusivity) increase with temperature. Fig. 
Q (a-b) shows the demagnetization timescales. We use 
an effective one-dimensional number of particles, which 
scales like where N is the total number of particles 
in a three dimensional trap. The fact that the demagne¬ 
tization rate due to the motion of the atoms in the trap 
can be comparable or to faster than the spin diffusion 
decay rate makes it clear that one relevant role of strong 
interactions is to suppress the effects of single-particle 
motion. 


IV. QUADRATIC MAGNETIC FIELD 

We now consider the case when the applied inhomoge¬ 
neous magnetic field varies quadratically with position: 
Bz{x) = mBx"^/A.fj,. Here B is a constant with units of 
frequency^. The total potential experienced by the atoms 
is 

V(x) = (w^ -I- (T^B) x^ . (18) 

T hus, each spin state sees a new trap frequency oja- = 
y/u}'^ + crB, with a the eigenvalue of . 

We can find the time dynamics of a harmonic oscillator 
eigenstate subjected to a sudden quench in the oscillator 
frequency by making the following ansatz, analogous to 
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Eq. § for the linear case: 


V’n {X,t) 


(19) 


where Oo- (t), jSa- (t), and (j)(y (t) real functions which are 
independent of n. It can be shown that this ansatz repre¬ 
sents a solution to the time-dependent Schrodinger equa¬ 
tion if the functions (t), Pa (t), and tpa (t) satisfy 



- 3 - 2-10 1 2 3 - 3 - 2-10 1 2 3 

x/bh x/a^ 


0.3 

0.2 

0.1 

0 

- 0.1 


a" (t) -I- -I- aB) aa (t) 


Pa (t) 

P'ait) 


eta (,t) ’ 

1 


cu 


(ty 


( 20 ) 


Equation (201 is a special case of the Ermakov equation, a 
well-studied nonlinear ordinary differential equation [Si- 
Defining B = BjuP as a dimensionless parameter, the 
solutions for our particular case are 


e^a (y) — 




Pa (t) = - 


2 -I- Ba -I- cos (2uit\/ 1 -I- Ba^ ^ 

2 (l -H Baj 

13^/1 + Basin (2ujt\/ 1 -I- BaJ 
2 -I- Ba -I- cos 1 -I- Ba^ ^ 


ujtVl + Ba 

Pa [t) = I - \Tr+ 


-I- — — arctan 


Vi + Ba cot (^ojt\/l + Ba)], (21) 


where [.J is the floor function. The solutions i/>„ (x, t) are 
wavefunctions whose probability densities maintain their 
original shape, up to periodic rescaling by a(t). The fre¬ 
quency associated with the “breathing” of the probability 
density is 2uj\/ 1 -I- Ba, and hence the periods of oscilla¬ 
tion are generally incommensurate for the spin up and 
spin down particles, see Fig. ([^. These modes are the 
breathing analogs of the well-known coherent states of a 
harmonic oscillator. 

Figure ([^ also shows 


{S^{x,t)) = i - |i/:;(x,t)p) , (22) 

that reflects the breathing of the spin up/down densities 
around the trap centers at x = 0. The different stan¬ 
dard deviations of the densities, also lead 

to a spin segregation, or spatial separation of the spin 
up/down densities: while the former becomes more lo¬ 
calized towards the center, the latter becomes more con¬ 
centrated at the edge of the trap. 


FIG. 5. Dynamics of spin observables in a quadratic 
magnetic field with B = O.BSca^. (a) {S^{x,t)) shows the 
spin up/down densities breathing in their traps centered at 
X = 0. The rates of their breathing, = \/u!^ ± B, are 
incommensurate. The standard deviations of the densities 
are hence spin up (down) becomes more con¬ 

centrated at the center (edge) of the trap, (b) Magnitude 
of (contrast), which measures the magnetization at 

position X. The magnetization decays when the spin up/down 
densities are separated. Note that decays/revivals of the mag¬ 
netization are aperiodic. 


A. Spin Dynamics 

To calculate (5+) = ^ J dx'tp^{x,t)ipi{x,t) we use 


Eq. (19) and write 


,A+^ _ i exp[z(^^(t) - pijt))] 
^ ^ 2 


(23) 


where 


L = 


[ dxe (&x) (cx) 


2”n 


a^{ty + a^{ty 


2a^(f)2a^(t)2 ' 


h = 




c = 




(24) 


and Hn{x) is a Hermite polynomial. Due to the difference 
in breathing frequencies, the spin up and spin down par¬ 
ticles accumulate a different dynamical phase, unlike the 
case of a linear magnetic field. Because of this dynamical 
phase difference, there is a net frequency shift. 

Using the generating function for the Hermite poly¬ 
nomials [ini, we can find a generating function for the 
integrals In as 


9i {z) = 


1 


•^/a — 2bcz + (6^ -I- — a) 




n—O 


(25) 


which leads to the following closed form for {S~^): 
(‘5'^) = d / T7Y(gig2)" X 


X 2 E 1 


2 ^aa^{t)ai{t) 

l-« -n (5?-l)(ffi-l) 

n ’ o ’ _2-2 


2 2 
9i9i 


(26) 
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Here, gi = h/y/a, g 2 = cl\fa^ and 2 F\ [.] is the hypergeo¬ 
metric function. To illuminate this result we can expand 
(.§+) to second order in ujt: 


(5+) 


2 L 


iB 


1 . 


1+73 (^ + n) + 


W 


^2 


(wt)^ -1-0 . (27) 


We see that at short times the signal acquires a frequency 


shift Ai/ = {B/uj){n+ and thus the collective spin ex¬ 
hibits a net precession in the X-Y plane. The frequency 
shift is a consequence of the additional frequency scale in 
the problem, namely, = y/uj'^ -\- B — — B Ki 

B/u) for B ^ Additionally, the contrast decays 
quadratically in Bt/uj, and depends on the mean of the 
squared energy. From the generating function, we can 
also find the Maxwell-Boltzmann thermal average exactly 
as gi (e“^“) /Z, which gives 


(^+) 


M.B. 


__ iexp ~ sinh (/3a;/2) __ 

y/(aj(t) + cosh(/3w) -f - /3;(t)) sinh(/3a;) - 2) 


(28) 


B. Spin-echo sequence for quadratic magnetic field 


To find the spin dynamics we solve Ermakov equation 
with the conditions of a spin-echo sequence. The result, 
expanded at short times, is 



• (29) 


As was the case for the linear magnetic field, the spin- 
echo removes the leading order contrast decay, and also 
removes the leading order frequency shift. We again in¬ 
terpret this result as due to the fact that the particles 
do not move to first order in time, and so the spin-echo 
can remove the effectively-static dephasing at lowest or¬ 
der. For longer times, there is an essential difference 
in the present case compared to the magnetic held with 
linear position dependence: here, {S'^) is not a peri¬ 
odic function of time, as the two frequencies + B 
and — B are incommensurate for general B. In the 
case of a linear magnetic held, the residual ehects of 
the spin-echo pulse on motion can be removed by ap¬ 
plying the spin-echo pulse after a single motional period, 
as both spin states oscillate at the trap frequency. In 
the quadratic case, the two spin states have different mo¬ 
tional periods, and so a spin-echo pulse cannot remove 
motional effects from both spin states at long times. 
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BTf 



FIG. 6. Quadratic magnetic field: (a) demagnetization time 
with and without spin-echo vs. temperature, where the Fermi 
temperature Tp = Nuj, and (b) contrast decay vs. time, 
using parameters from Ref. m- The timescales are on the 
order of or faster than the timescale of spin segregation in the 
experiment (~ 100ms). Additionally, application of a spin- 
echo pulse leads to faster demagnetization. 


C. General Considerations 


For the quadratic magnetic held the contrast decay 
due to motional ehects, expressed as a demagnetization 
rate, also exhibits a quadratic and quartic scaling with 
time with and without a spin-echo pulse, respectively. 
However, unlike the linear magnetic held case, there are 
generally no full revivals of the magnetization, due to 


the incommensurate breathing frequencies of the spin up 
and spin down densities. Using the approximation of a 
continuum of energies and a constant density of states, 
the corresponding demagnetization timescales are given 


















by: 


field case at zero temperature (Ref. [S]). 


^ ^ ) (E2 U?ji) 

TT^ ln(e^^‘^ - 1) + ln(e^^‘^ - 1)^ - dLig ( y =^) 
3/33w(iV +1) ^ 

(30) 

where the last line gives the mean squared energy. In 
comparison to the case of a constant gradient, the rate is 
proportional to the mean squared energy rather than the 
mean energy, so the demagnetization rates increase more 
quickly with particle number and temperature than the 
linear case. In Fig. we plot the demagnetization time 
and contrast decay for the quadratic magnetic field case 
with and without spin-echo, using experimental parame¬ 
ters from Ref. El , but with the proper rescaling of parti¬ 
cle number to account for the different dimensionality, as 
was done for the linear magnetic field case. Under these 
conditions a spin-echo pulse again leads to faster demag¬ 
netization. We find that the demagnetization timescales 
are on the order of or faster than the spin segregation 
timescales observed in an experiment in which interac¬ 
tions were non-negligible (~ 100ms). There, the segre¬ 
gation timescales were dictated by the mean interaction 
energy of the gas. The faster short-time decay exhib¬ 
ited by the non-interacting system reveals once more the 
suppression of motional effects from interactions. 
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V. SUMMARY AND OUTLOOK 

We have demonstrated that the motion of non¬ 
interacting fermions in a trap can lead to demagneti¬ 
zation on timescales faster than those caused by in¬ 
teractions in recent experiments. Additionally, a spin- 
echo sequence can have the counterintuitive effect of en¬ 
hancing the rate of demagnetization at times compara¬ 
ble to or longer than the motional period. The anal¬ 
ysis present here, which exactly characterizes all the 
relevant timescales and parameters that determine the 
non-interacting spin dynamics of a finite temperature 
fermionic gas, can serve as a platform for a better un¬ 
derstanding of the interplay between motional-induced 
and interaction-induced spin dynamics. 

A logical direction for future work is to characterize 
the crossover from non-interacting to interacting dynam¬ 
ics. A promising avenue along these lines is to treat 
the occupation of single-particle energy states as fixed 
and consider the effect of interactions on these occupied 
single-particle states. This will be explored in detail in 
an upcoming work. 

During the finalization of this work we became aware 
of a complementary paper that treats the linear magnetic 
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Appendix A: Derivation of generalized 
coherent-state formnla 


Consider H = ^ + = uj a + ^) and suppose 

that 14)) is a solution to the (time-dependent) Schrodinger 
equation. We now show that j^/^) := D (z := zoe“*“‘) If/)) 
is also a solution, where D is the phase-space displace¬ 
ment operator D{w) = exp (wa4 — w*a) and zq is any 
complex number, corresponding to the initial displace¬ 
ment in phase space. 

Observing that [a,D{z)] = zD{z) and [a''', D{z)] = 
z*D{z), we have 


[H, D{z)\ = oj ^za^£>(z) -I- z*D{z)a^ 

= OJ (za^ + z*a — z*z) D{z). 


(Al) 


Noting also that z = —iujz,z* = iojz*, we have from 
the Baker-Campbell-Hausdorff formula. 



dt\ J 

-lOJza'^Diz) + e^“^(-fwz*a)e-^*“e-^o^“/2 
(—zwza'^ — ioJZ*a + iujz*z) D{z), 

(A2) 


since [e^“\d] = —ze^“\ Thus 

= w {za) + z*a - z*z) D{z)\(j)) + D{z) 

= [iL,i)(z)]|<(>)+.D(z)iL|<(.) 

= HD{zM = H\ij). 

(A3) 







